The two dimensional J\f = (2, 2) Wess-Zumino Model in the Functional 
Renormalization Group Approach 
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We study the supersymmetric N = (2, 2) Wess-Zumino model in two dimensions with the func- 
tional renormalization group. At leading order in the supercovariant derivative expansion we recover 
the nonrenormalization theorem which states that the superpotential has no running couplings. Be- 
yond leading order the renormalization of the bare mass is caused by a momentum-dependent wave 
function renormalization. To deal with the partial differential equations we have developed a nu- 
merical toolbox called FlowPy. For weak couplings the quantum corrections to the bare mass found 
in lattice simulations are reproduced with high accuracy. But in the regime with intermediate cou- 
plings higher-order-operators that are not constrained by the nonrenormalization theorem yield the 
dominating contribution to the renormalized mass. 
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I. INTRODUCTION 

In the search for high energy theories beyond the stan- 
dard model supersymmetric models are a topic of great 
interest. Supersymmetry reduces the hierarchy and the 
fine-tuning problem. It has to be broken at some energy 
scale since supersymmetry has not been observed in low 
energy physics. The breaking does not occur on the per- 
turbative level and therefore nonperturbative tools are 
needed to analyze these models, e. g. lattice formulations 
or the function renormalization group. 

Lattice formulations and simulations have been suc- 
cessfully applied to nonperturbative problems in field 
theory. Although there has been considerable progress in 
the last years [11-13 there are still some difficulties in the 
lattice formulation of supersymmetry. The discretization 
leads to a (partial) breaking of supersymmetry and the 
implementation of dynamical fermions on the lattice still 
poses a challenge. 

Nonperturbative continuum methods, such as the func- 
tional renormalization group (FRG) which manifestly 
preserve supersymmetry, can complement the lattice cal- 
culations. The FRG has previously been applied to a 
wide range of nonperturbative problems such as criti- 
cal phenomena, fermionic systems, gauge theories and 
quantum gravity, see e. g. [5l4l3| for reviews. Applied to 
supersymmetric theories it circumvents problems of the 
lattice formulation such as supersymmetry breaking due 
to discretization. But in order to solve the FRG equa- 
tions truncations have to be employed which introduce a 
different kind of error. 

Quite a few conceptual studies of supersymmetric the- 
ories in the framework of the FRG have already been per- 
formed. The main ingredient is the construction and use 
of a manifestly supersymmetric regularization scheme. 
For example such a regulator has been presented for the 
four-dimensional Wess-Zumino model in [14, Hal. Inves- 



tigations for one-, two- and three-dimensional M = 1 
Wess-Zumino models have been performed in [l6l - [l9j . 
A FRG formulation of supersymmetric Yang-Mills the- 
ory employing the superfield formalism has been given 
in 20] ; for further applications see also [2l|, [22[ . General 
theories of a scalar superfield including the Wess-Zumino 
model were studied with a Polchinski-type RG equation 
in which yields a new approach to supersymmet- 
ric nonrenormalization theorems. The nonrenormaliza- 
tion theorem has also been proven with FRG methods in 
[23l | . In [24j a Wilsonian effective action for the Wess- 
Zumino model by perturbatively iterating the FRG is 
constructed. 

The aim of this work is twofold. On the one hand, 
we want to compare the results from the supersymmetric 
formulation of the FRG equations to lattice data for the 
renormalized mass in the two-dimensional Af = (2, 2) 
Wess-Zumino model [I). This comparison allows us to 
estimate the truncation error. The renormalized mass is 
defined as the location of the pole of the propagator in the 
complex plane therefore we have to take the momentum 
dependence in the FRG framework into account. 

There are several applications, where this dependence 
is important but the related contributions lead to a much 
higher numerical effort for the solution of the flow equa- 
tions. Full momentum dependence of propagators and 
vertices has previously been treated successfully in the 
literature [25l432j |. We introduce a numerical toolbox 
called FlowPy which is designed to solve the flow equa- 
tions with momentum dependence as encountered in this 
paper. FlowPy can be adapted to solve not only flow 
equations with momentum dependence but also other 
differential equations encountered in the FRG framework 
e. g. for field dependent effective potentials. In this paper 
we demonstrate that FlowPy solves the flow equations 
reliably. 

The paper is organized as follows: In section [TT] we 
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introduce the Af — (2, 2) Wess-Zumino model in two 
dimensions. In section IIIII we sketch the derivation of 
the supersymmetric flow equations for the superpotential 
and the (momentum-dependent) wave function renormal- 
ization. The flow equation for the superpotential will 
lead to the nonrenormalization theorem. In section IIVI 
first FlowPy is described and then we specialize our flow 
equations and demonstrate that perturbation theory is 
reproduced correctly. In section[V]we compare the renor- 
malized mass calculated in the FRG approach with the 
results from lattice simulations. 



II. THE N = (2, 2) WESS-ZUMINO MODEL IN 
TWO DIMENSIONS 



The system has two bosonic ground states which lead to 
a nonzero Witten index (35j . therefore supersymmetry 
is never spontaneously broken in the Af — (2, 2) Wess- 
Zumino model. 

A characteristic feature of the Af = I Wess-Zumino 
model in four dimensions survives the dimensional reduc- 
tion, namely that bosonic and fermionic loop corrections 
cancel in such a way that the effective superpotential 
receives no quantum corrections. This is called the non- 
renormalization theorem [36 - 38] . In the two-dimensional 
model considered here the cancellations even render the 
model finite. The model has been studied intensively in 
the literature, see e. g. [1, 0, [39l - l4]| for lattice simula- 
tions. 



The Af = (2, 2) Wess-Zumino model in two dimensions 
can be found by a dimensional reduction of the Af = 1 
model in four dimensions [33j. The Lagrange density is 
given by 

Jzf = 2d4>d(t) + 4>Mil) - -FF + -FW'{(j>) + -F W'(4>) 

2 2 ^ 

with Dirac fermions ip and -0. The fermion matrix M 
reads 



M = $ + w"(4>)p+ + w"{4f)P- 



(2) 



with P± = (1 ± 7*)/2, F = F l + iF 2 and = 0!+ icf) 2 
as well as d = \{d\ — id 2 ) and z = X\ + ix 2 . The super- 
potential is denoted by W(4>) — u{<f>i,<j> 2 ) + w{4>i,4>2)- 
We work in the Weyl basis with 7 1 = 01, 7 2 = — 02 and 
7* = i7 X 7 2 = 03. The complex spinors can be decom- 
posed as ip = (ipi 1P2) and if> = Oipx ^2) ■ The Lagrange 
density is invariant under the supersymmetry transfor- 
mations 



8<f> = ipiSi + sitpi, S(f> = ip 2 e 2 + £2^2, 

Stpi = —r^Fei + d(j)e 2 , Si/; 2 = dfoi - ^Fe 2 , 

SF = 2{&4> 1 e 2 - e 2 dipi), SF = 2{dfp 2 ei - sidtp 2 ). 



(3) 



The superspace formulation of this model is constructed 
in appendix [X] A detailed discussion of the underlying 
supersymmetry algebra and a construction of the super- 
space can be found e. g. in |34| . 

Integrating out the auxiliary fields yields the on-shell 
Lagrangian 

C on = 2d4>dc/) + hv'{<t>)W'{<t>) + tpMi/j. (4) 
In this paper we will consider the superpotential 



W{ct>) = \n<? + \gf. 



(5) 



III. SUPERSYMMETRIC RG FLOW 

Following the lines of our previous works [l6l - [l9| we 
construct a manifestly supersymmetric flow equation in 
the off-shell formulation. Our approach is based on the 
FRG formulated in terms of a flow equation for the effec- 
tive average action T^, i.e. the Wetterich equation [42| 



d k V k = - STr 



.(2) 



Rk 



dkRi 



(6) 



The scale dependent T k interpolates between the micro- 
scopic action S for k — > A, with A denoting the micro- 
scopic UV scale, and the full quantum effective action 
r = ■ As the model considered in this paper is UV- 
finite, the cutoff A can be set to infinity. The interpolat- 
ing scale k denotes an infrared (IR) regulator scale below 
which all fluctuations with momenta smaller than k are 
suppressed. For k — > 0, all fluctuations are taken into ac- 
count and we arrive at the full solution of the quantum 
theory in terms of the effective action T. The Wetterich 
equation defines an RG trajectory in the space of action 
functionals with the classical action S serving as initial 
condition. 

The second functional derivative of T k in Eq. © is 
defined as 



fr (2) "| = JLr k — 

V k Jab d^ a 0* 6 



(7) 



where the indices a, b summarize field components, in- 
ternal and Lorentz indices, as well as spacetime or mo- 
mentum coordinates. In the present case, we have 
\E fT = (4>,<f>, F, F,ip,ip) where * is not a superfield, but 
merely a collection of fields. The momentum-dependent 
regulator function R k in Eq. ([6]) establishes the IR 
suppression of modes below k. In the general case, 
three properties of the regulator Rk(p) are essential: (i) 
Rk(p)\ p 2/k 2 ^o > which implements the IR regulariza- 
tion, (ii) Rk(p)\k 2 /p 2 ^o = which guarantees that the 
regulator vanishes for k — > 0, (iii) Rk{p)\k^A->-oo — > 00 
which serves to fix the theory at the classical action in 
the UV. Different functional forms of Rk correspond to 
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different RG trajectories manifesting the RG scheme de- 
pendence but the end point T k ^ — > T remains invariant, 
see e. g. Refs. [1, 143l - l48j |. Supersymmetry is preserved 
if the regulator contribution to the cutoff action ASk 
(cf. Eq. (|lip ) is supersymmetric. 

As an ansatz for the effective action we use an expan- 
sion in superspace (see Appendix [Al for our conventions j3 



r fc = 



d 2 x 



dydy Zl(dd)^ 



(8) 



2 / d 2 x / dy W fe ($) - 2 / d 2 x / dy W k {$) 



4^ 



Z^p z )[2p z ^ + ^0~-FF 



(9) 



In contrast to the usual supercovariant derivative ex- 
pansion we have included those combinations of the su- 
percovariant derivatives that merely reduce to spacetime 
derivatives. A momentum dependence in W k is irrelevant 
as found in section UlIBI An arbitrary Kahler potential 
$) integrated over the whole superspace) is not 
taken into account here, since we expect only a small 
influence for the renormalized mass. Another contribu- 
tion neglected in this truncation comes from the terms of 
higher than quadratic order in the auxiliary field and the 
corresponding supersymmetric partner terms, denoted as 
auxiliary field potential. In the following we will only 
work with real and imaginary part <f>i, <f>2, Fi, F%. 

For this scale dependent effective action the auxiliary 
fields obey the equations of motion F — W k (4>) / Z 2 and 
F = W' k ((j))/Zl. This leads to the on-shell action 



1 k — 



4tt 2 



1 

+ 2 



Z 2 k (p 2 ) 



-iZ 2 k ( P 2 )^ + HW k l P + 



(10) 



A. Supersymmetric regulator 



Supersymmetry is preserved if we shift the mass by 
a momentum-dependent infrared regulatoi0, m — > m + 
Z 2 -ri(k,p 2 ) or multiply the wave function renormal- 
ization by a momentum-dependent regulator function, 
Z 2 — > Z k - r 2 (fc,p 2 ). Such regulators are the same as 
the ones used in the previous models EMU. To get a 



1 For the Fourier transformation we use the convention dj 



with the notations p = (pi,P2) T and p = \p\ where there is no 
risk of misunderstanding. 
2 The regulator function is multiplied with the wave function 
renormalization to ensure reparametrization invariance of the 
flow equation. 



regularized path integral R k is included in terms of the 
cutoff action AS k . It reads in a matrix notation 



AS k = l 



— - * zIrI ^ t 



(11) 



with * = (</>i fa Fx F 2 ip{-p) T ip(j>)) and 




!'h ( 'n '' ,','/ with R% = 



and Rf. 



p 2 r 2 -1 n ■ cr 3 

n • <t 3 -r 2 ■ 1 
ip ■ r"2 — T\ ■ 1 

r 2 + n ■ i o 



(12) 



With these regulators at hand we can proceed to calcu- 
late the flow equation. Inserting ansatz ^ in the flow 
equation ([6]) , the propagator can be calculated along the 
lines described in [49| : The fluctuation matrix rjj, 2 ' + R k 

(2) 

is decomposed into the propagator Ig + R k including 
the regulator functions and a part Ar^. containing all 
field dependencies. The flow equation ^ is expanded in 
the number of fields, see Appendix [B] for the expansion 
and the explicit matrices. 



B. Flow equation for the superpotential 
nonrenormalization theorem 



The 



The quantity at leading order in the supercovariant 
derivative expansion is the scale dependent superpoten- 
tial. We obtain the flow equation by projecting onto the 
terms linear in the auxiliary fields. We can choose either 
the real or imaginary part of the auxiliary field as they 
are bound to give the same results due to supersymme- 
try. The superpotential W{4>) = u{4>\,4>2) + iv(<j>i,<j>2) 
is a holomorphic function of <j>\ and cj>2, and therefore 
its real and imaginary part obey the Cauchy-Riemann 
differential equations 



du 

Mi 



dv 



du 



Mi 



Using these equations we find for the flow equations of 
the superpotential 



d k u k = 0, d k v k = =* d k W k = d k W k = 0. 



(14) 



This means that the superpotential remains unchanged 
during the RG flow. The Kahler potential does therefore 
not influence the flow of the superpotential, as found in 
[3?| . Even the nontrivial momentum dependence consid- 
ered here does not change this result. The nonrenormal- 
ization theorem is verified by the flow equations in the 
present truncation. This result is similar to the proofs in 
four dimensions discussed in [231 ] and [111 ]. 

As the flow vanishes at leading order, the first quan- 
tity with a nonvanishing flow is the wave function renor- 
malization which is a term at next-to-leading order in 
our truncation. It will turn out later that the momen- 
tum dependence is important for the renormalized mass 
(cf. Sec. El therefore we already include it in ansatz ©. 
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C. Momentum-dependent flow equation for the 
wave function renormalization 

The flow equation for the wave function renormaliza- 
tion can be obtained from a projection onto the terms 
quadratic in the auxiliary fields. It is derived in Ap- 
pendix |B] and reads 



d k zl{p) = -t 



d 2 q h{p~q)h(q) 
Aw 2 v(q) 2 v(p — q) 



[dfciJi (q - p)M(p - q)v(q) + d fe i?! (q)M (q)v{p 
2 f d 2 q h(p- q)d k R 2 (q)u(q)v(p - q) 



q)\ 



•V 



4tt 2 



v(q) 2 v(p — q) 2 



d 2 q h(q)d k R 2 (q - p)v(q)u(p - q) 



Am 2 v(q) 2 v(p — q) 2 

with the abbreviations (recall that |q| = q) 

Kq)={r 2 (q) + l)Z 2 (q), 
M(q) =m + n(q)Z 2 (q) , R % {q) = n (q) Z 2 (q) , 
u(q) =M(q) 2 - q 2 h 2 (q), v(q) = M{qf + q 2 h\q) 



(15) 



(16) 



Here we are dealing with a UV-finite theory and there- 
fore it is sufficient to use the simple, masslike infrared 
regulator 



r\(k,p 2 ) — k and r 2 (k,p 2 )=0. 



(17) 



After a shift in the integration variables in the second 
part of the integral (fl"5|) the flow equation simplifies to 



d 2 q kZ\ (q) + m 



d k Z 2 (p) = -16g 2 J 



At: 2 N(q) 2 N(p-q) 
Z 2 (q)Z 2 (\p-q\)d k (kZ 2 (q)), (18) 

where we have introduced the abbreviation 

N(q) = (q 2 Zt(q) + (kZ 2 ( q ) + m) 2 ). (19) 

In order to deal with the partial differential equa- 
tion we have developed a (parallelizable) numerical tool- 
box called FlowPy. In the next section we present 
the main ideas behind our numerical setup to solve the 
momentum-dependent flow equation. A detailedpresen- 
tation of FlowPy is deferred to a future paper [50| • 

IV. NUMERICAL SETUP 

Structurally, the flow equation to be solved numerically 
is of the form 



d k z{k;p) = J d n ql[k]p]q;z{k;fi{p,q))] 

z(k;f 2 {p,q));...;z(k;f n (p,q))] , (20) 

where the external momentum p is treated with a dis- 
cretized grid. While the integrand / may actually also 




FIG. 1. Perturbative flow for the parameters A = g/m — 0.3 
and m = 1 . The solid line is the plot of Eq. (|2ip . 



be a function of d k z(k;p), and hence the integral flow 
equation be given in implicit form only, numerical re- 
sults suggest that, at least in the model studied here, 
d k z(k;p) should be sufficiently small to allow an itera- 
tive approach, where the integrand is evaluated first un- 
der the assumption d k z(k;p) = 0, and the result is then 
used to re-evaluate the integrand with a better approxi- 
mation to d k z(k;p) until convergence is reached. Apart 
from this conceptual issue, the technology to deal with 
an evolution equation of this kind is readily available in 
an accessible form via the SciPy "Scientific Python" ex- 
tension [5l| to the Python [52[ programming language. 
The details of the numerical strategy are described in 
Appendix [Cj 

As a test for the numerical approximation and the abil- 
ities of FlowPy, we solve the flow of the perturbative wave 
function renormalization. It is inferred by setting Z k (q) 
to its classical value Z k (q) = 1 on the right hand side 
of Eq. (fT8|) . The perturbative flow with 60 discretization 
points is shown in Fig. Q] for different values of k. 

It is possible to calculate the perturbative expres- 
sion for Zf_ loop (p) analytically from the perturba- 
tive flow equation by performing the fc-integral using 
lim/j^oo 7*i, r 2 — > oo and linu^o r i, r 2 — > 0. This yields 



7 2 

-^l-loop 



9 



d 2 q 



(m 2 + q 2 )(m 2 + \q — p\ ) 



, artanh (p(Ar 



Trp^/Ar 



- p A 



(21) 



which is shown as a solid line in Fig. [T] This shows 
that possible errors in the numerical calculation of the 
wave function renormalization with FlowPy due to dis- 
cretization, interpolation and the boundary condition 
Z k (q —> oo) = 1 are under control. We will consider Z k 
and the renormalized masses obtained from it as exact in 
the employed truncation. 

In the next section we will determine renormalized 
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masses from the nonperturbative wave function renor- 
malization with full momentum dependence calculated 
with FlowPy. 



THE RENORMALIZED MASS 



The analytic continuation of the bosonic propagator 



Gbos(j?) 



1 



P 



T,(p,m,g) 



(22) 



has a pole which defines the renormalized mass. Since 
the bare mass to is a parameter of the superpotential ([5]) 
it is not changed during the flow. E is the self-energy. 
As expected from a supersymmetric theory, the pole of 
the fermionic propagator leads to the same renormalized 
mass as the bosonic propagator. 

The Fourier transformation of G DO s(p) yields the cor- 
relator 



dp 
2^ 



G{ Pl ,oy pixi . 



(23) 



The renormalized mass can be obtained from the long 
range exponential decay of this quantity and is in the 
following denoted as correlator mass TO corr . One can also 
define a renormalized mass, which we denote as propaga- 
tor mass, through m 2 Top = (G bo s(p)) _1 |p=o • 

To compare the renormalized masses from the FRG 
with the results of the lattice simulation we consider 
the masses of the particles in the on-shell theory. In 
our truncation the bosonic propagator from the on-shell 
action (fTU)) reads in the infrared limit 



Gbos (p) 



1 



and the fermionic propagator reads 



Gfcrm (P) 



P 2 Z£^ (P 2 ) + m 2 



(24) 



(25) 



Both propagators have the same poles and therefore lead 
to the same renormalized masses for bosons and fermions. 

For the propagator mass the fields in the on-shell ac- 
tion have to be rescaled with the wave function renor- 
malization such that the kinetic term is of the canonical 
form. Neglecting the momentum dependence in the wave 
function renormalization we obtain 



''prop 



(26) 



For a small self-energy E a comparison between Eq. 
and (|24| leads to the approximate relation 



zl_ 



•o(P) 



£(p,ra,g) 
p z — 77 r 



(27) 



A numerical calculation can provide Z 2 (p) only for 
real p and its analytic continuation cannot be deter- 
mined straightforwardly. Instead we consider the dis- 
crete Fourier transformation of G^ (p) with momenta 
p = {0, 2n/aN, 2tt(A - l)/aN} on the interval x G 
[Q,aN = L\, For distances much smaller than L this 
should approximate C^ s ° (x) in a well-defined way. More 
precisely, instead of the exponential decay one gets the 
long distance behavior 



C a , mcor (a;i) oc cosh(m CO rr(a;i - L/2)) 



(28) 



after the integration over the spatial direction. The mass 
can be determined from a fit to this function, as it is done 
in lattice simulations. The details of this procedure can 
be found in Appendix [D] 

With the analytic result (|2"Tj) for Z 2 _ loop at hand we 
can calculate the poles of G^®{p) and obtain a per- 
turbative approximation of m corr .Note that this analytic 
solution of the perturbative flow together with eq. (|27|) 
leads to the same result as a one-loop on-shell calcula- 
tion of the polarization E (cf. Appendix |E|) . Expanding 
the pole of the propagator (|22|) to first order in the di- 
mensionless parameter A 2 = g 2 /m 2 leads to the one- loop 
approximation of the renormalized mass 



(to 



c Vr p ) 2 =- 2 (l-^A 2 + 0(A 4 )) . (29) 



A. Weak couplings 

Let us start with an investigation of the weak coupling 
sector which is defined as A < 0.3, where perturbation 
theory provides an excellent cross-check to establish the 
correctness of our ansatz and the errors in the numerical 
approximation. 

The bare mass in the lattice simulations Q is taken 
to be to = 15. Concerning the units of the mass note 
the following: In the lattice calculation, the mass is mea- 
sured in units of the box size, i. e. the physical volume of 
the lattice simulation. Similarly, everything can be refor- 
mulated in terms of the dimensionless ratio of bare and 
renormalized mass. 

For the numerical treatment of Eq. (|T5|) we have to 
use dimensionless quantities. Because of the nonrenor- 
malization theorem the bare mass quantities in the su- 
perpotential enter in the flow equation only as parame- 
ters. Rescaling the dimensionful quantities with the bare 
mass sets the scale in this model. We have set this scale 
to to = 1. To get the same units as in the lattice simula- 
tions the resulting renormalized mass is multiplied with 
15. 

The correlator masses in the weak coupling regime are 
calculated with the momentum-dependent wave function 
renormalization from the flow equation (|18l) solved with 
FlowPy. The technical details of the determination of 
the correlator masses are described in Appendix [D] The 
results are shown in the second column of Table HI The 
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14.4 



0.05 0.1 0.15 0.2 0.25 0.3 

g/m 

FIG. 2. Comparison between lattice data taken from [3] and 
our results for the correlator mass with momentum 

dependence and mj^p without momentum dependence in the 
weak coupling regime 



TABLE I. Renormalized masses obtained with the flow equa- 
tion with and without momentum dependence (m™? an d 
m prop) as weu as lattice data m^" 100 from a continuum ex- 
trapolation [3] in the weak coupling regime. 



A 


m FRG 


FRG 
/ < f p r op 


lattice 
" ^corr 


0.02 


14.998 


14.998 


14.999(1) 


0.04 


14.991 


14.992 


14.993(3) 


0.06 


14.979 


14.983 


14.977(4) 


0.08 


14.963 


14.970 


14.963(5) 


0.10 


14.943 


14.952 


14.935(6) 


0.12 


14.917 


14.931 


14.905(9) 


0.14 


14.888 


14.907 


14.871(9) 


0.16 


14.854 


14.878 


14.83(1) 


0.18 


14.815 


14.846 


14.82(1) 


0.20 


14.773 


14.810 


14.75(2) 


0.22 


14.674 


14.771 


14.71(2) 


0.24 


14.674 


14.728 


14.63(2) 


0.26 


14.619 


14.681 


14.60(2) 


0.28 


14.559 


14.631 


14.53(2) 


0.30 


14.496 


14.578 


14.45(3) 



values in the fourth column are taken from a Monte-Carlo 
simulation on the lattice |]| . We discuss the lattice results 
further in Sec. IV Bl For the time being it suffices to note 
the agreement of lattice and perturbative results within 
the statistical errors. Hence perturbation theory already 
provides a good check for our results. 

In Fig. [2] we show the correlator masses from the flow 
equation, the lattice simulation and the one-loop result 
(f2"9"| for m corr . The masses calculated from the flow 
equation agree very well with perturbation theory and 
with the results from lattice simulations. This can be 
quantified by comparing the correction to the bare mass 
Am CO rr = m- m corr . We find Am™/Am^ icc ~ 0.95. 



Taking into account the statistical error of the lattice 
data no significant difference to the FRG results can be 
found. 

We conclude that in the weak coupling regime the trun- 
cation of the flow equation with full momentum depen- 
dence suffices to capture the main aspects of the model. 
Higher-order-operators, which yield an auxiliary field ef- 
fective potential, have, as expected, little influence. 

To investigate the influence of the momentum depen- 
dence in the wave function renormalization, we calculate 
the propagator mass (l26|) . The results are shown in the 
third column of Tab. U and in Fig. [2] A comparison be- 
tween the propagator mass and the correlator mass from 
the lattice calculation yields Am™/Am[ a ' r f c ~ 0.75. 
Already in the weak coupling regime it is necessary to 
include the momentum dependence in order to determine 
the corrections to the renormalized mass with satisfying 
accuracy. 

B. Intermediate couplings 

At intermediate couplings 0.3 < A < 1 a significant de- 
viation of our numerical results from perturbation the- 
ory can be observed. In this regime the perturbative 
calculations can no longer provide a reliable test for the 
numerical results and we have to rely on lattice calcu- 
lations. In a supersymmetric theory their result must, 
however, be considered with care as a lattice formulation 
of supersymmetry still poses difficulties [53[ . A common 
approach is to implement only a part of the supersym- 
metry which allows to recover the complete symmetry in 
the continuum limit in many cases. 

In the present model there are further complications 
for a lattice formulation, especially in the intermediate 
coupling sector [4|. The considered discretizations are 
invariant under half of the supersymmetry They suffer, 
however, from the dominance of a contribution to the 
action that is a mere discretization of a surface term at 
larger couplings. In general the correct continuum limit 
can only be obtained with unrealistically high numeri- 
cal effort. The relevance of this effect depends on the 
coupling strength and on the specific discretization. For 
intermediate couplings the nonlocal SLA C discretization 
and the Twisted Wilson discretization provides the most 
reliable results (cf. [4] for details). The renormalized 
masses of these discretizations are used for a comparison 
with our results. They are shown in the third and fourth 
column of Table HTlFI and displayed in Fig. [3] (boxes with 
error bars) together with the order A 2 expanded result 
(f2"9")l for m corr (dashed line). 

Note that the spontaneous breaking of the Z2 sym- 
metry introduces also finite volume effects in the lattice 
simulations. Although there are well known prescriptions 



3 All lattice results are extrapolated to the continuum. 
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FIG. 3. Comparison between lattice data taken from [3] and 
our results for the correlator mass with momentum 

dependence and without momentum dependence in the 

intermediate coupling regime 



TABLE II. Masses obtained with the flow equation with and 
without momentum dependence m™? and m FR p) as well as 
lattice data [J] in the regime with intermediate couplings. 



A 


m FRG 

' ' "corr 


m FRG 

' ' "prop 


tw. Wilson 


SLAC imp. 


0.35 


14.321 


14.428 


14.23(2) 




0.40 


14.123 


14.259 


13.99(3) 


14.00(1) 


0.45 


13.905 


14.069 


13.62(5) 




0.50 


13.666 


13.861 


13.30(6) 




0.55 


13.411 


13.636 


12.8(1) 




0.60 


13.138 


13.394 


12.2(1) 


12.44(6) 


0.65 


12.854 


13.137 


11.9(2) 




0.70 


12.556 


12.866 


10.4(5) 




0.75 


12.248 


12.583 






0.80 


11.932 


12.290 




10.2(3) 


0.85 


11.609 


11.987 






0.90 


11.280 


11.676 






0.95 


10.948 


11.358 






1.00 


10.613 


11.036 




8-1(311 



C. Wozar, private communication 



to implement these properties of the theory in lattice sim- 
ulations, they still lead to additional complications |54j . 

The correlator masses determined from the FRG are 
shown in the second column of Tab. [TT] and displayed 
in Fig. [3] (crosses). Additionally the perturbative result 
for the renormalized mass is shown (solid line), which is 
determined from the pole of the propagator (|2"2"j) with the 
polarization calculated in Appendix [En . 

Although the corrections from the wave function renor- 



malization with full momentum dependence to the bare 
mass capture some of the quantum effects, they do not 
account for all the nonperturbative effects present in this 
model. To quantify this, we compare these corrections to 
the corrections found in lattice calculations. This yields 
results between Am™/A 



and ATO™P/Am' atticc 



m™™ ce ~ 0.9 for A = 0.35 
,, , ;n _,,.,.,„.,. _ 0.65 for A = 1.0. The fact that 
the wave function renormalization accounts for less of 
the quantum corrections as the coupling grows is due to 
the growing influence of higher-order operators especially 
the auxiliary field potential. In the present truncation we 
have only considered terms that are at most quadratic in 
the auxiliary field and have neglected back reactions from 
a potential for the auxiliary field. As can be seen from a 
diagrammatic expansion of the flow equation, terms up to 
order directly modify the flow equation for the wave 
function renormalization, which is proportional to Ff . It 
is known from our previous investigations of scalar super- 
symmetric models fl6l | that the influence of higher order 
operators grows with the strength of the couplings. A 
truncation that goes beyond the momentum-dependent 
wave function renormalization has to be considered to 
improve the results in the regime with intermediate cou- 
plings. 

The results for the propagator mass are shown in the 
third column of Tab. Ql] and in Fig. |3] (triangles). Com- 
pared to the lattice results we find Am^p/Am^'' 106 ~ 
0.75 for A = 0.35 and Am^/Am^™ ~ 0.6 for 
A = 1.0. The improvement due to the momentum de- 
pendence in Z\ is not as pronounced as it is in the weak 
coupling regime. 



VI. CONCLUSIONS 

In this paper we have applied the functional renormal- 
ization group to the M = (2, 2) Wess-Zumino model in 
two dimensions. The model is UV-finite which allows a 
direct comparison to results from lattice simulation. 

The first quantity to be calculated in a supcrcovari- 
ant derivative expansion is the superpotential. It is well 
known from the nonrenormalization theorem that it does 
not receive quantum corrections. In the language of the 
FRG the nonrenormalization theorem is recovered in a 
very simple form, namely that the superpotential has a 
vanishing flow equation. The proof only uses the fact 
that the superpotential is a holomorphic function and 
therefore the Cauchy-Riemann differential equations for 
its real and imaginary parts hold. 



4 The perturbative results in this regime have to be interpreted 
with care. The dashed line is an expansion of the on-shell one- 
loop pole mass to order 0(A 2 ) whereas the solid line is the result 



of an off-shell calculation. Both results agree up to order A 2 
but perturbation theory can no longer be trusted in the regime 
with intermediate coupling strength. The on-shell calculation 
has to fail at large values of A 2 because otherwise the renormal- 
ized masses will become negative. To preserve supersymmetry 
in the RG flow the FRG uses an off-shell formulation. Therefore 
it is not unexpected that it is close to the off-shell perturbation 
theory. 
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Hence the first term in the expansion that receives 
a correction from renormalization is the wave function 
renormalization. It leads to the renormalization of the 
bare mass in the on-shell theory, with the renormalizcd 
mass defined as the pole of the propagator in the complex 
plane. We have calculated the renormalized mass with 
and without momentum dependence in the wave function 
renormalization. 

In order to benchmark our results we use lattice cal- 
culations. In the weak coupling regime the results for 
the renormalized mass calculated without the full mo- 
mentum dependence capture only 75% of the quantum 
corrections to the bare mass whereas 95% of the correc- 
tions are captured if the full momentum dependence in 
the wave function renormalization is taken into account. 
This leads to the conclusion that the momentum depen- 
dence of wave function renormalization dominates in this 
regime. Higher-order-operators only have a small influ- 
ence. 

For intermediate couplings the picture changes. We 
have investigated the complete contribution to the flow 
from the momentum-dependent wave function renormal- 
ization. Our findings are that in this truncation only 
65% of the quantum corrections to the bare mass de- 
termined in the lattice simulations are captured for the 
largest coupling considered in this paper. Without mo- 
mentum dependence 60% of the corrections are gener- 
ated. This leads to the conclusion that in the regime 
with intermediate couplings the momentum dependence 
in the wave function renormalization does not include 
all important contributions to the renormalized mass. 
Instead, the quantum corrections generated by higher- 
order-operators which lead to an auxiliary field potential 
are expected to be relevant for the renormalized mass. 
They have to be included in order to reduce the devi- 
ations between the results from lattice calculations and 
the FRG. The calculation of these contributions as well 
as contributions from the Kahler potential remains an 
interesting challenge for future work. 

Although there have been great improvements in the 
simulations of the model on the lattice they still suffer 
from finite size effects and the finite lattice spacing, which 
leads to a breaking of supersymmetrjjf]. An interesting 
application of the FRG is an analysis that includes finite 
volume effects which could help to estimate the influence 
of the finite size. This can allow one to separate it from 
the discretization errors. 

The analysis presented in this paper can easily be ap- 
plied to the Af = 1 Wess-Zumino model in four dimen- 
sions from which the two-dimensional Af = (2, 2) model 
is derived. Especially the nonrenormalization theorem 
for the superpotential emerges in the same way. In both 



models the momentum-dependent wave function renor- 
malization is the first relevant contribution in the covari- 
ant derivative expansion and the flow equations differ 
only in the measure of integration. As we have found in 
the present model the effective potential for the auxiliary 
field is expected to dominate the quantum effects as the 
strength of the coupling constant grows. 

For the treatment of the partial differential equation 
we have developed FlowPy, a numerical toolbox. It can 
be applied in quite generic situations to solve the FRG 
equations and to calculate contributions such as the full 
momentum dependence of vertices. We hope that nu- 
merical software like FlowPy will help to obtain better 
predictions from FRG calculations. We would like to 
add that it can also be applied for the calculation of an 
arbitrary potential V(<p) instead of the full momentum 
dependence Z(p). We plan to make FlowPy available 
soon 1501. 
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Appendix A: Superspace formulation 

The superspace formulation is constructed from the su- 
persymmetry transformations. The chiral and antichiral 
superfields can be obtained from the lowest component 
by acting on it with the exponentiated supersymmetry 
transformations 

4 i 

$0, z, a, a) = exp(-6 a )<f>(z, z) = V — {-5 a ) n (f>(z, z) 

=4>(u, u) — -01 (u, u)ai — aiipi(u, u) — — a±a± 

4 f 

$(z, z, a, a) = exp(-8 a )(/>(z, z) = J~] —J-5 a ) n (f)(z, z) 

L — ' TV. 

n=0 

- _ F(u, u) _ 
=4>(u, u) - ijj 2 (u, u)a 2 - ai2ip2(u, u) ^ — a 2 a 2 

(Al) 

with the chiral variables u = z — ^a 2 a\ and u = z + 
\ol\ol 2 . The supercharges are 



5 In contrast to common lattice formulations the basic require- 
ments of locality and reflection positivity are broken in the cur- 
rent simulations of this theory. This was done to reduce the 
unavoidable breaking of supersymmetry on the lattice. 



n 9 , 1 * n 9 X - « 

oa\ 2 oai 2 

O 9 ^ 1 A h 9 l ~ ft 



da 2 2 



da 2 2 



(A2) 
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and the supercovariant derivatives read 

£>i = -tt= ^«20, Di = h -a 2 o 

doti 2 dai 2 . . 

5 1 9 1 - ( j 

D 2 = -— L> 2 = -=. h-ai9. 

oa2 2 aa2 2 

The superfield obeys the (anti)chiral constraint 

D 2 $ = D 2 <£> = 0, Di$ = L>ii> = 0. (A4) 
The supersymmetry transformations are generated by 

5$ = (eQ + Qe)$, «J$ = (eQ + Qe)#. (A5) 
The Lagrange density is given by 

C =£kin + £pot (A6) 

= - 2 /" dydy$$-2 J dy W($) - 2 J dy W(§) 
with dy = (iaidai and (iy = da-ida.2- 



Appendix B: Flow equation for the 
momentum-dependent wave function 
renormalization 

To obtain the flow equations for the wave function 
renormalization we decompose the second derivative of 
the effective action into a field independent part Iq +R k 

(2) 

and a field dependent part AT k (in the following we 
drop the momentum dependence of the regulators for 
simplicity of notation): 

^(2) 



(T^+R k )(q,q')+AT k (q,q') 



A 



S(q-q') 



AA AC 
AD AB, 



with (h ■ 
A 

and 



(l + r 2 )Z 2 (q),M=(r 1 Z k 2 (q)+m)) 

Ml 



q 2 h-t M-a 3 
M-a 3 -h-t 



Bn = H. 



(Bl) 



(B2) 



AA = 2g 



AC = 2g 



AD 



AB = 



2g 
2.9 



/ Fx -F 2 


0i 


-02 \ 


-F 2 -Fx 


-02 


-01 


01 -02 








\-02 -01 





o / 








i>i -#2 













\o o J 






( ipi #i 


o o\ 


(<? + q 




o oy 




'01 + «02 







I « 


/>i - i 



(B3) 



The flow equation can then be expanded [49j in 
d t T k = ift STr ((if) + TJkJ^AT) 

-ia t STr((rf + i? fc )- 1 Ar 



(B4) 



with dt acting only on the regulator. STr denotes a trace 
in field space as well as an integration in momentum 
space. The wave function renormalization is a term pro- 
portional to F 2 and can be obtained from the second 
term in this expansion. To calculate this we define 

M(q,q')= f (4 2) +R k )- 1 (q)6(q + q")Ar(q",q>) 

Jq" 

=(T^+R k )- 1 (q)AT(-q,q') (B5) 
and the second term in the expansion reads 

9 t Str/ M(q,q')M(q',q) 
•> <?,<?' 

= Str [ (if > + R k )-\q)d t R k {Tf ) + R^iq) 
J q,q' 

x AT(-q,q')(T^ +R k )- 1 (q')Ar(-q',q) 
+ Str f (T^+R k )- 1 (q)AT(-q,q')(T^+R k )- 1 (q') 



x d t R k (q')(T^ +R k )- 1 (q')AT(-q',q) (B6) 

where Str denotes a trace in field space. We take the 
functional derivative with respect to F(p) an d F(— p) 
and set all fields to zero in order to project on the wave 
function renormalization Z k (p 2 ). This yields 



d k Z 2 k {p) 



2 f d2 1 h(p-q)h(q) 



4ir 2 v(q) 2 v(p — q) 2 
d k Rx(q-p)M(p-q)v(q) + d k R 1 (q)M(q)v(p-q)} 



d 2 q h(p - q)d k R 2 (q)u(q)v(p - q) 



4tt 2 



v(q) 2 v(p — q) 2 



4g 
Ag 

Air 2 v(q) 2 v(p — q) 2 

with the abbreviations 



, , d 2 qh(q)d k R2(q-p)v(q)u(p-q) 



(B7) 



Kq) = (ra (q) + 1) Z 2 (q) , 
M(q) =m + n(q)Z 2 (q) , R % (q) = n («) Z 2 (q) , (B8) 
tt(g) =M(q) 2 - q 2 h 2 (q) 1 v(q) = M(qf + q 2 h 2 {q) 

Appendix C: The computational strategy 

To solve flow equations depending on an external mo- 
mentum the computational strategy is as follows: 
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• Mapping the k —> flow to a forward time evolution 
problem by introducing k = —k. 

• Discretization of the problem by approximating z(k;p) 
with an interpolation function that is determined by N 
support points. These are best chosen to be equidistant 
on a logarithmic scale, accounting for the expectation 
that a wide range of scales should contribute in a com- 
parable way to the integral. 

• For now, we have been using interpolation in con- 
junction with the two-dimensional integration function 
scipy. integrate .dblquadO from Scientific Python. 
However, this ad-hoc approach should allow great ef- 
ficiency improvements by instead considering a tighter 
integration of numerical quadrature with adaptive dis- 
cretization. In particular, it should then also be pos- 
sible to also pass on information about the discretized 
flow equation's Jacobian matrix to the numerical ODE 
integrator. There hence is considerable potential for 
further efficiency improvements of numerical RG flow 
code. 

• Solving the resulting ordinary differential equation 
for the values of z{k;p) at the support points with 
SciPy's scipy . integrate . odeint () function (which 
internally uses lsoda from DDEPACK). 

Concerning the numerical solution of the ODE, some 
manual tweaking of integration parameters such as max- 
imal step sizes is required when the coupling constant g 
is large and k ~ m. We note that, by the very na- 
ture of this problem, the computation is readily paral- 
lelized: The effort to numerically determine the right 
hand side integral is expected to roughly grow like 0(N 2 ) 
with the number TV of support points, and computations 
for different support points are independent. In com- 
parison to the computational effort required to compute 
the integrals, the communication overhead to distribute 
the values of z(k;p) at different support points is fairly 
negligible, hence using one of the readily available MPI- 
extensions to Python to intelligently distribute the work- 
load becomes an attractive option. One should, however, 
take care that the core structure of the integrand then 
is implemented in a compiled (C code) Python extension 
before even thinking about parallelization. 



Appendix D: Determination of the renormalized 
mass 



The numerical calculations of Z 2 in the main text use 
a grid of N = 60 points in the direction of p 2 , dis- 
tributed equidistantly on a logarithmic scale. The result 
for Zk^,o{p) is interpolated with splines to calculate the 
propagator G^°(p). A discreet Fourier transformation 
of Gbos°(P) yields the correlator C(xi) on the interval 
x\ e [0, L] with n = 10 001 intermediate points. In the 
main text we have used L — 15. From its large distance 



behavior 



C a ,m aoi ( x i) « cosh(m corr (xi - L/2)) 



(Dl) 



the correlator mass m corr is determined by a least 
square fit. The fit range is constrained to the interval 
[%i,skip, ■ ■ ■ , L — x\ t ski P \ where the contributions of ex- 
cited states are negligible. The value of Xi :S ki P is deter- 
mined such that Tn coa (xi g kip) shows a plateau. We can 
either fit on the whole range [x\ yS kip, — xi, s kip] ~ 

this quantity is called mf°^ &l - or make the fit just in- 
side a small interval of size 0.2 starting from cci jS kip _ this 
quantity is called m l °£f. 

In the left panel of Fig. @] m]?^ 1 is shown for two dif- 
ferent discretisations of Z\{p 2 ), N — 200 in the upper 
and N = 600 in the lower panel. In the right panel the 
same is shown for mfj,°^ al . From these plots we can read 
off that for x\ yS kip not too large there is a clear plateau 
which is stable if the box size is increased. But for very 
large Xi lS ki P the local mass oscillates. As this oscillation 
is reduced if the discretization is increased it is due to 
fluctuations in the spline interpolation of Z^.. At small 
values of the correlator the numerical errors are more 
important for the masses. As the fluctuations become 
visible for large box sizes, in these cases the global mass 
fit is of no use because it averages over the local mass 
and is strongly influenced by the oscillations. We will 
therefore take the plateau of rn^f as the value of the 
renormalized mass. 



Appendix E: Perturbation Theory 

In perturbation theory the mass is determined from 
the one particle irreducible (proper) vertex The 
perturbative expansion of this vertex is 



o + o 



(El) 



0(g 4 



e(p) -2. 9 2 y 



2 f d 2 s 6m 2 — 4((p — s) 2 + to 2 ) 



2^ 



47T 2 (s 2 + TO 2 )((s — p) 2 + TO 2 ) 

d 2 s 4(s 2 -s-p- ie^s^{s v - p v )) 
Air 2 (s 2 + m 2 )((s - p) 2 + to 2 ) 



/ d 2 s 2(to 2 -p 2 ^ 



4-7T 2 (s 2 + m 2 )((s — p) 2 + TO 2 ) 



=2g 2 / dz 
Jo 

=2g 2 [ dz- 



2(m 2 ~p z ) 



At: 2 (s 2 + z(l — z)p 2 + m 2 ) 2 
2(to 2 -p 2 ) 



2\2 



j 47r(z(l — z)p 2 + TO 2 ) 
4g 2 (m 2 — p 2 ) 



~Kp\J 4m 2 + p 2 



artanh (p(4m 2 + p 2 ) _1/2 ) (E2) 
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FIG. 4. Left panel: m l °™ with A = 0.6 for discretizations N = 200 and N = 600 and different box sizes L — 15, 25, 35 and 45. 
Right panel: with A = 0.6 for discretizations TV = 200 and TV = 600 and different box sizes L = 15, 25, 35 and 45. 
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